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Abstract. An analysis of Hq and H/3 spectra in a sample of 30 cool dwarf and subgiant stars is presented using 
MARCS model atmospheres based on the most recent calculations of the line opacities. A detailed quantitative 
comparison of the solar flux spectra with model spectra shows that Balmer line profile shapes, and therefore the 
temperature structure in the line formation region, are best represented under the mixing length theory by any 
combination of a low mixing-length parameter a and a low convective structure parameter y. A slightly lower 
effective temperature is obtained for the sun than the accepted value, which we attribute to errors in models 
and line opacities. The programme stars span temperatures from 4800 to 7100 K and include a small number of 
population II stars. Effective temperatures have been derived using a quantitative fitting method with a detailed 
error analysis. Our temperatures find good agreement with those from the Infrared Flux Method (IRFM) near solar 
metallicity but show differences at low metallicity where the two available IRFM determinations themselves are i n 
disagreement. Comparison with recent temperature determinations using Balmer lines by Fuhrmann (1 



1998 



2000) 



who employed a different description of the wi ng absorption due to self-broadening, does not show the large 
differences predicted by Barklem et al. (2000b). In fact, perhaps fortuitously, reasonable agreement is found 
near solar metallicity, while we find significantly cooler temperatures for low metallicity stars of around solar 
temperature. 
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>■ 1. Introduction 



Hydrogen is by far the most abundant species in typi- 
cal stellar atmospheres. In late- type atmospheres, a small 
fraction of hydrogen atoms capture free electrons form- 
ing ions, which dominate the continuum opacity. 
Hydrogen itself is the main opacity source for earlier spec- 
tral types. The few lines that its simple atomic structure 
makes in the spectrum have a very distinct sensitivity to 
the atmospheric properties compared to metal lines. In op- 
tical stellar spectra, absorption lines of the Balmer (n=2) 
series are commonly used to study photospheres. The well 
populated lowest levels of the atom produce considerable 
opacity at the centre of the lines, and interactions with 
charged ions, electrons and other hydrogen atoms result 
in extended wings in high-density atmospheres. In late- 
type dwarfs, these wings are believed to form very close to 
LTE, in the deepest photospheric layers. As most protons 
are bound to electrons forming hydrogen, and hydrogen 



influences the main continuum opacity source, changes in 
the hydrogen abundance are hardly reflected in the lines' 
strengths, and the strengths are much more weakly af- 
fected by gravity and metal abundances than perturba- 
tions to the temperature. These properties attracted the 
attention of stellar spectroscopists, making Balmer lines a 
key feature in stellar classification schemes (e.g. Morgan 
et al. |1943|). 



Detailed analyses of 
stars are more recent (e 
al 



et 



in late-type 
Fuhrmann 

1994 |1994|; van't Veer-Menneret et al 



Balmer lines 
g. Gehren 



1981 



1998 



Gardiner et al. 1999). These modern studies exploit 



progress in theory and experiment on line broaden- 
ing to infer stellar effective temperatures from Balmer 
lines. Vidal et al. ( |1970| , |1973D developed a success- 
ful unified theory to model the interaction of hydrogen 
atoms with charged particles. Those calculations have 
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been recently superseded by Stehle (1994) and Stehle 
& Hutcheon ( 1999D , who have computed Stark broad- 
ened line profiles including ion dynamic effects under the 
model microfield method. A further important broaden- 
ing contributor is the collis ions w ith neighbouring hydro- 
gen atoms. Ali & Griem (196£) used the multipole ex- 
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pansion of the resonance interaction potential in the im- 
pact approximation to calculate line-widths from this pro- 



cess. Barklem et al. (2000a, 2000b, hereafter paper I and 
II respectively) have presented a self-broadening theory 
accounting also for dispersive-inductive interactions and 
without use of the multipole expansion. It was shown that 
the new description of the self-broadening would have a 
large impact on the computed Balmer line profiles, in par- 
ticular when applied to derive effective temperatures for 
metal-poor dwarf stars. 

Assuming a proper understanding of the line broaden- 
ing of the hydrogen lines, the most notable difficulty for 
the use of their wings as a temperature indicator is the 
fact that they are formed in very deep layers. The ther- 
mal structure of the deepest optically transparent layers 
in late-type stars is significantly affected by convection. 
Simple modelling of surface convection is still a challenge. 
The commonly used mixing-length formalism incorporates 
unphysical parameters which are hard to connect with 
quantities that can be derived from observations or hy- 
drodynamical simulations and, therefore, are difficult to 
constrain. This obstacle makes flux-constant homogeneous 
models particularly uncertain in these layers. Other theo- 
ries that dispense with the free parameters in the mixing- 
length theory (MLT) have also been proposed (e.g. Canuto 
et al. 1996|; Canuto & Mazzitelli 1991); however, Gardiner 
et al. ( 1999| ) found that after adjusting the mixing-length 
parameter a, MLT performed simila rly. S emi-empirical 
modelling (e.g. AUende Prieto et al. 2001 ) does not of- 
fer a viable solution, as the employed metal lines do not 
probe layers as deep as those where hydrogen lines are 
formed. These theoretical problems related to establish- 
ing a model atmosphere combine with observational con- 
straints. Balmer lines require high dispersion observations 
with a large spectral coverage, and a predictable instru- 
mental response, making possible a methodical and ac- 
curate continuum normalisation. In this situation, use of 
Balmer lines as a part of spectroscopic analyses requires 
a careful assessment of all possible sources of error. 

In this paper we investigate the impact of the afore- 
mentioned line broadening theory advances in the frame- 
work of ID model atmospheres, and attempt to iden- 
tify those areas where improvement is most desirable. In 
Sect. 2 we describe our observations and reduction proce- 
dure. In Sect. 3 we describe how model spectra are com- 
puted, and put the atmospheric models in context with 
others. In Sect. 4 we introduce a method for quantita- 
tive comparison of observations and model spectra, and 
subsequently an automated fitting procedure for deriving 
effective temperatures. In Sect. 5 we make a detailed sur- 
vey of possible errors and their effect on effective tem- 
peratures. Application to the solar spectrum and the pro- 
gramme stars is then presented. Finally in Sect. 6 the re- 
sults are compared with other work, and in Sect. 7 our 
conclusions are presented. 



2. Observations and Reduction 

The majority of the targets were observed at the 
2.5 m Isaac Newton Telescope (INT) at La Palma 
with the MUSICOS cross -dispe rsed echelle spectrograph 
(e.g. Baudrand & Bohm 1992) fed by a fibre from the 
Cassegrain focus. Targets were selected in order to obtain 
a sample of dwarf and subgiant stars with a spread of 
temperature and metallicity, with emphasis on well stud- 
ied stars. Due to the constraints of the telescope used, only 
a relatively small number of metal-poor targets could be 
observed. The data consist of observations taken in May 
1999 and January 2000, using 2048x2048 and 1024x1024 
pixel CCD arrays respectively. The spectra have a resolu- 
tion of approximately R = X/6X w 30000, and signal-to- 
noise ratio (SNR) per pixel of typically better than 100 at 
Ha and H/3, reaching 300 for some of the brighter targets. 
The Ha and H/3 spectra have an average SNR of 180 and 
150 respectively. The 1999 observations were obtained as 
a backup programme during periods of high cirrus. We do 
not expect the subsequent scattering to be significant and 
no attempt to remove the water vapour lines is made. 

With this spectrograph, the Balmer lines are well cen- 
tred in the orders; however, the broad wings of these lines 
can span several orders, making determination of the con- 
tinuum level a difficult process. This may be resolved by 
merging several consecutive orders, but a good normalisa- 
tion is needed. Flat-fielding of MUSICOS spectra taken at 
the INT has long been a problem, as the internal flat-field 
lamps available at the Cassegrain focus were designed for 
instruments with low and medium spectral resolution. The 
MUSICOS spectrograph is a high-resolution instrument, 
and the front end of the fibre is relatively small compared 
to the entrance slits of the lower resolution instruments, 
thus missing a significant portion of the light from the 
lamps. The flat-field frames obtained during our observa- 
tions were not satisfactory, seen for example by the fact 
that the behaviour of the continuum of orders when di- 
vided by the flat-fleld frames required a high-order poly- 
nomial to fit them, rather than a low-order polynomial 
typical of the residual difference from the stellar flux dis- 
tribution and the flat-flcld lamp. 

For broad lines the pixel-to-pixel variations in the CCD 
response function are not as important as for narrow lines 
and so instead of using traditional flat-fields, we con- 
structed artiflcial flat-flelds which neglect these variations 
but reconstruct the blaze shape from stellar exposures. 
The relation between the blaze shapes of the different or- 
ders is in principle a smoothly changing function, and be- 
cause in general only a few orders contain broad spectral 
lines the blaze shape for these orders can be determined 
by interpolation. The spectral layout of frames taken with 
the MUSICOS spectrograph is very well adapted to this 
approach, because not only is the number of orders in the 
frame large (of order 50) , but also the coverage in the dis- 
persion direction is large (typically 60-lOOA), tracing the 
shape of the blaze function far out in its wings. Orders 
have significant overlap in wavelength, even at Ha where 
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the overlap is smallest, the amount depending on the CCD 
array used. To determine the shape of the blaze function 
we used the fact that for our data set the extracted pro- 
files contain only absorption lines. We then performed a 
polynomial fit to the upper points in each extracted order. 
This method fails if there are no continuum points avail- 
able due to the presence of a broad hydrogen line in the 
order. A surface fit of the neighbouring orders is then used 
to determine the artificial flat-field at any problematic or- 
ders. The regions for interpolation are selected interac- 
tively, which may include other problematic regions such 
as strong metal lines or molecular bands. This procedure 
can also be used if the frames are first divided by a tra- 
ditional flat-field exposure, as the fiat-field blaze-function 
is still smoothly changing. On the larger CCD fringing 
at Ha made this procedure necessary for the rectification 
of this line. However, for other cases we did not divide 
the exposures as we found the variation of the orders was 
in general smoother. Orders are extracted and scattered 
light corrected using standard IRAF echelle package rou- 
tines, before normalisation as described above using IDL. 
Finally the wavelength calibration is determined, orders 
are merged together weighted by the SNR at the wave- 
length in the order, and wavelengths are corrected to the 
the stellar rest frame. The method is illustrated in Fig. ^ 
for one of the more difficult cases (low SNR«120 and 
broad line). 

For Procyon, in addition to spectra from 
INT/MUSICOS, a very high quality spec trum has 



been obtained by AUende Prieto et al. (2002) with the 
McDonald Observatory 2.7 m telescope and the 2dcoude 



spectrograph (TuU et al. 1995) for a large spectral region 
including H/3. The spectra have SNR of about 1000 
at H/3, and a resolution of approximately 2 x 10^ (see 



AUende Prieto et al. 2002 for details). The same contin- 
uum placement techniques described for INT/MUSICOS 
were applied to these observations. 

The spectrum of HD 103095 (Gmb 1830) was acquired 
on 28 and 29 April 1990 with the McDonald 2.7 m tele- 
scope and the Coude spectrograph. The observations em- 
ployed a conventional grating in first order and a colour 
filter to eliminate light from higher orders. The detec- 
tor was a 800x800 pixel CCD. The dispersion is about 
0.12 A/pixel and the SNR was a few hundred. This spec- 
trum was processed using IRAF being normalised using a 
second-order polynomial fit. 

Comparison of reduced spectra from different expo- 
sures of the same targets, revealed similar results to those 
of Fuhrmann et al. ( 1993| ), namely an internal consistency 
of approximately 1% of the continuum flux, dependent on 
SNR. Most of this error probably stems from the deter- 
mination of the continuum placement. The error can be 
reduced through co-adding a number of independent ob- 
servations (Fuhrmann et al. 1993| ). Direct compari son o f 
moon spectra with the FTS atlas (Kurucz et al. 1 



Comparison of the McDonald and INT Procyon spectra 
shows similar agreement, despite quite different instru- 
mental setups. We will test this consistency later by ex- 
amining scatter in temperatures derived from different ob- 
servations. 

3. Models and Synthetic Spectra 

In this work we employ ID LTE plane-parallel mod el at 
mospheres from the MARCS code (Asplund et al 



1997 



version, see this reference and those therein for details). 
The adopted version of MARCS employs opacity sam- 
pling techniques, except in the infrared where opacity 
distribution functions are deemed sufficiently accurate. 
This circumvents problems of opacity distribution func- 
tions namely allowing free c hoice of individual metal abun- 
dances (cf. Fuhrmann et al. 1997 ). Convection is modelled 
under MLT which is parameterised in terms of a number of 
parameters, the most important of which for model atmo- 
spheres and our discussions are the mixing- length I, which 
is usually expressed in units of the pressure scale height 
Hp as the mixing-length parameter a = I /Hp, and the 
structure parameter y describing the temperature struc- 



ture within convective elements (see Henyey et al. 1965 
for details). A logTioss (tioss is Rosseland mean optical 
depth) grid suitable for computation of Balmer lines was 
chosen. The depth grid extends in to logTioss — 2 and the 
grid is finest in the region where these lines are formed 

-1 < logTross < 0.6. 

Models are computed with scaled-solar photospheric 
abundances from Anders & Grevesse (1989) except for 
and most importantly the iron abundance where a "low" 
value of log(A^Fc/A^H) + 12 = 7.50 is adopted in agreement 
with the meteoritic value (e.g. Grevesse & Sauval 1998). 
A microturbulence of 1.5 km s~^ was used in all cases 
for computing the models. The MLT parameters used 
were a = 1.5, y = 0.076 the default MARCS s etting , 
and a = 0.5, y = 0.5 following Fuhrmann et al. ( 1993|) . 
Hereafter we will use MARCS to refer to the former and 
MARCS05 to the latter. The reasons for these choices are 
discussed later. Each grid was also recomputed with en- 
hancement of the alpha-elements by 0.4 dex, which we will 
denote in discussions with an "a" (e.g. MARCS05a). When 
referring specifically to the grids with scaled-solar abun- 
dances we will use an "s" . The model grids are computed 
with spacings of 100 K in T^s, 0.1 dex in log g and 0.5 dex 
in [Fe/H], within which we interpolate specific models. 

Synthetic model flux spectra are computed as de- 
scribed in paper II. In short, spectra are computed as- 



suming LTE using SYNTH (Piskunov[l992|), Stark broad- 
ening is described by the model-microfield method calcu- 
lations of Stehle & Hutchcon ( 1999D and self-broadening 
calculations of paper II are used. Radiative broadening is 
included, as is an estimate of the helium coUisional broad- 



degraded to appropriate resolution shows agreement at 
the level of 0.5-1%, although small differences in wave- 
length calibration make the direct comparison difficult. 



ening (based on rescaling of Barklem & O'Mara 1997), al- 
though both effects are small. We will refer to this as the 
STEHLE-I-BPO broadening recipe. For comparison with 
previous work we also do calculations substituting the 
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Fig. 1. An example of the continuum determination procedure, for the H/3 line in HR 5447, obtained on 29 May 1999 
with INT and MUSICOS. The upper panel shows the extracted order which is centred on H/3. The smooth lines show 
the fit to the continuum for nearby orders where full lines are determined from the spectra, and dashed lines were 
determined by interpolation from the full lines. The thicker dashed line is the continuum corresponding to the plotted 
spectral order. Note the pixel number is arbitrary as the ends of orders have been removed to avoid poor behaviour 
of the polynomial fits. The rapid variation of the blaze function across an order is a characteristic of the Littrow 
configuration used in MUSICOS. The lower panels show the final spectrum (with best fit model spectrum) along with 
an estimate of the SNR per pixel in the continuum which also serves to show the position of each order. The upwardly 
convex regions are the centres of the orders, and the concave parts the overlapping ends of the orders which have 
higher formal SNR due to contributions from the two overlapping orders. For clarity the region corresponding to the 
order centred on the Balmer line (i.e. that spectrum shown in the top panel) is shaded, with the areas of overlap with 
the adjacent orders shaded darker. Note the cosmic ray hit in the upper panel has been removed. 



Ali & Griem (1966) resonance broadening theory for the 
self-broadening calculations, everything else remaining the 
same. This will be referred to as the STEHLE+AG broad- 
ening recipe. In our implementation of the Ali & Griem 
theory we only consider the broadening of the 2p state, 
since the 2p-nd transition dominates (see discussion in pa- 
per II). Note that this differs from Fuhrmann et al. (1993) 
where both the 2p and np state are included. 



3.1. Comparison to other models 

To put the MARCS05 models which will be used in most of 
this work in context, we now compare them differentially 
with some other models. To compare the models particu- 
larly as regards computed Balmer line profiles, we made 
comparisons by finding the Tiff's of the MARCS05 Balmer 
line profiles which best match profiles produced by other 
models using least-squares minimisation. This should in- 
dicate approximately the difference in Toff that would be 
found from the input model relative to MARCS05 for the 
given line. STEHLE+BPO recipe was used to produce all 
profiles and only the wings are considered. 
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Table 1. Differences between other solar models, in 
terms of effective temperatures from Balmer lines, of our 
MARCS05 models. 



Model 



GS 



(Grevesse & Sauval |1999| ) 
HM 



(Holweger & Miiller |197J) 
MISS 



(AUende Prieto et al. gOOl|) 
KOVER 

(Kurucz 



1993, see paper II) 



KNOVER 



(K) 


(K) 


-81 


-196 


-114 


-192 


-132 


-250 


-24 


-214 


+73 


-10 



(Kurucz 1993, see paper II) 



First we compared with other common solar models, 
both semi-empirical and Kurucz (1993) theoretical flux- 
constant models. Table || shows the difference in temper- 
ature of the MARCS05 model needed to match the in- 
put profile (i.e. Toff(best match)— 5777) for the first two 
Balmer lines. Large differences with semi-empirical models 
are seen, which produce much weaker Balmer lines than 
the MARCS05 solar model. When comparing with solar 
observations the semi-empirical models produce lines too 
weak and the MARCS05 too strong when STEHLE-fBPO 



is employed (see paper II and Sect. 5.2). The differences 
with Kurucz models in Table ^ stem from different MLT 
parameters and inclusion of overshoot in the KOVER 
model. We should point out that after publication of paper 
II it was realised that the MARCS model used in that pa- 
per was in fact computed with a = 1.25 and y — 0.5, the 
same as the Kurucz models, not a = 1.5 and y — 3/(47r^) 
as stated, hence the reasonable agreement of the pro- 
files with KNOVER. One should note that despite the 
rough agreement of H/3 temperatures from KNOVER and 
MARCS05, the line shapes are markedly different, which 
was true of H/3 in all cases. 

The differences with the Munich group models 



(see Fuhrmann et al. [1997| , hereafter MUNICH) are of 
interest to understand the discrepancies between the 
Tcff values from this work and Fuhrmann (1998, 2000| ). 
Comparisons were made for a range of stellar temper- 
atures and metallicities and the results are shown in 
Table ^. Alpha-enhancement is included for all models be- 
low solar metallicity. A single variation of gravity is also 
shown. We see that the effect on Ha is generally larger 
and can be as high as 40 K for hotter stars. For most cases 
the corrections are positive, which indicates a hotter tem- 
perature would be found using MARCS05 than MUNICH 
models. Fig. ^compares temperature structures from both 
model grids. We see that in the upper regions of the atmo- 
spheres, MARCS05 models typically have slightly hotter 
and less steep T-t structures. We note that for the cooler 
models, particularly at low metallicity, the temperature is 
hotter throughout the entire atmosphere. 



Table 2. Differences between MARCS05 and MUNICH 
models in terms of effective temperatures derived from 
Balmer lines. 



Model 


Areff(Ha) 


Areff(H/3) 


roff/log3/[Fe/H] 


(K) 


(K) 


5200/4.2/0.0 


+6 


+3 


5600/4.2/0.0 


+9 


+11 


6000/4.2/0.0 


+21 


+13 


6400/4.2/0.0 


+39 


+1 


5200/4.6/0.0 


+3 


-2 


5600/4.6/0.0 


+6 


+ 1 


6000/4.6/0.0 


+16 


-1 


6400/4.6/0.0 


+29 


-6 


5200/4.2/-1.0 


-8 


-15 


5600/4.2/-1.0 


-1 


-11 


6000/4.2/-1.0 


+ 14 


-7 


6400/4.2/-1.0 


+28 


-6 


5200/4.2/-2.0 


+2 


-8 


5600/4.2/-2.0 


+2 


+1 


6000/4.2/-2.0 


+ 11 


+1 


6400/4.2/-2.0 


+31 


+2 







; -.-^^^ [Fe/H] 


= 0, log g = 4.2 




50 

















-50 








40 


[Fe/H] 


= -2, log g = 4.2 2 




20 






h- 











-20 







-3 -2-10 1 

log ^5000 

Fig. 2. Differences between T-t structures of MARCS05 
and MUNICH models, for T^s = 5200, 5600, 6000, and 
6400 K (full, dotted, dashed and dot-dashed lines respec- 
tively), all for log 5 — 4.2. The upper and lower panels 
show [Fe/H]= 0.0 and [Fe/H]= —2.0 cases respectively. 



4. Fitting Method 

Quantifying the comparison between a synthetic profile 
and an observed profile, allows us to automate the fitting 
process through minimisation of the chosen statistic, and 
also provides a statistical measure of the goodness-of-fit. 
This further allows the subjectivity associated with such 
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fitting, usually done by inspection, to be shifted as we will 
demonstrate. We employ a reduced statistic, namely 



1 



X 



N~ M 



N 

E 



(1) 



where N is the number of wavelength points, Af is the 
number of free parameters (here one, namely T^s), fi is the 
synthetic residual flux, Fi the observed residual flux, and 
CTj = 1/SNR. For the value of a-i = 1/SNR we estimated a 
constant average value of the SNR (in the continuum) for 
the given Balmer line observation. Despite the fact that 
the formal SNR of our observations varies across the line as 
shown in Fig. |^, the regions with highest SNR are those at 
the end of the orders where continuum placement may be 
less certain. Therefore these should not be given a higher 
weight, and a constant average value is more appropriate. 

For our case of Balmer lines, blending lines must be 
accounted for. Calculations including the metal lines are 
impractical as current spectral line databases have not 
yet reached a level of completeness or accuracy where all 
observed lines in the solar spectrum in the Balmer line 
regions are well accounted for (see Fig. 9 in paper II). 
Furthermore, such an approach would require an iterative 
procedure with more free parameters. 

A much simpler approach is to attempt the identifi- 
cation of spectral windows across the line which are ex- 
pected to be free from blends in the stars of interest. This 
was done by using the high resolution FTS solar spectrum 



(Kurucz et al. 1984) and the best of our spectra of hot- 
ter and cooler stars of around solar metallicity (namely 
Procyon and HR 8832) to identify windows which would 
be free from obvious blending across most of our sample. 
High resolution spectra of a metal-rich K dwarf would be 
particularly useful in this respect, as this perhaps repre- 
sents the hardest case, but to our knowledge such spectra 
are not available. We make the assumption that unblended 
spectral windows do in fact exist and that we see them in 
the above mentioned high quality spectra. This is reason- 
able for the sun and Procyon at least at Ha and H/? based 
on examination of the spectra, but may for other stellar 
parameters or higher Balmer lines be questionable. Once 
these windows had been identified we considered other 
desirable qualities of the windows as a set. Firstly the 
Balmer line cores, which are formed high in the atmo- 
sphere in non-LTE are excluded, only the pressure broad- 
ened wings are considered. We attempted to give approx- 
imately equal weight to all parts of the line wings, despite 
a natural tendency for the inner-wings to show less blends 
due to saturation by the hydrogen line. We note that if the 
windows are evenly distributed across the line wings, and 
a least-squares fitting is employed such as in this work, 
the procedure is quite similar to matching the equivalent 
width of the wings of the line, and the information on 
the line shape is not biased towards any part of the line. 
Furthermore, a reasonable number of windows need to be 
employed so as to avoid the possibility that statistical fluc- 
tuations in the observations affect the fltting procedure. 



One should also be sure that windows are appropriate for 
the spectral resolution of the observations, and this will 
be discussed in more detail shortly. 

When comparing model spectra with observations in 
this way, we employ the same spectral windows, or a 
slightly reduced subset, for all stars. Rejection of some 
windows is necessary in specific cases, such as the fast- 
rotating early F stars where metal lines are very broad 
and encroach on these unblended regions, or if an atmo- 
spheric line by chance falls within a window in the refer- 
ence frame of a given star, or if a window happens to fall in 
the line core. This method restricts the human interaction 
in the fitting process, once the initial windows are decided, 
to simply this choice of rejecting some windows. By fol- 
lowing this approach we attempt to avoid the situation 
where one must distinguish between noise fluctuations in 
the spectra and possible weak metal lines, which can be 
particularly difficult in metal-poor stars with lower SNR 
observations. Thus in essence, the subjectivity involved 
in the usual fitting by visual inspection is shifted to this 
choice of windows. 

In cool stars, where self-broadening is important, one 
must be cognizant that both descriptions of this process 
used in this work employ the impact approximation. The 
extreme limit of validity for the self-broadening in the 
STEHLE+BPO recipe is given by 



AAmax = 3.654 X lO^^A^ 



(2) 



where the result AAmax is in A, with A the central wave- 
length in A, T the temperature in K, aiQi the cross-section 
in atomic units for v = lO^m s^^ and a the dimensionless 
velocity parameter (not to be confused with the MLT a 
parameter). Typical values for the sun are discussed in 
paper II. The region of validity varies through the stellar 
atmosphere, but based on the assumption that the line 
wings are predominantly formed in the region of the at- 
mosphere around optical depth unity, a reasonable esti- 
mate for a given stellar profile can be obtained by set- 
ting T = Toff - Spectral windows falling outside this region 
of validity are rejected. For consistency, even when the 
STEHLE+AG recipe is used we use the same limit as for 
STEHLE +BPO , although a calculation based on the AH 
& Griem ( |1966| ) cross-section would give a slightly larger 
region of validity. 

If errors in the data are normally distributed, minimi- 
sation of the statistic provides the maximum likelihood 
estimate of the model parameters, in our case T^g. Our 
errors will in fact be Poissonian, approaching a normal 
distribution at high SNR. A more robust (insensitive to 
departures from normal distribution) statistic to better 
account for outliers might be of benefit. 

For each star we computed a grid of Hq and H/3 pro- 
files for a given model grid at 10 K intervals in Tcn around 
an initial guess based on literature values. The statis- 
tic was always found to be well behaved varying smoothly 
with Toff, approximately parabolic in form with only a 
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Fig. 3. The fitting method for solar Ha and H/3 profiles, here showing the best fit found with MARCS05 and 
STEHLE+BPO which corresponds to the parameters in Table |5| The shaded regions show the windows used for 
determining the statistic. The full vertical lines show the estimated limit of validity of the impact approximation. 
Note the windows outside this region are rejected and thus no residual is plotted. 





Fig. 4. As in Fig. || for Procyon Ha and H/3 profiles. 
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single minimum (no other local minima). We therefore 
searched this grid for the best fitting profile, and finally 
performed a polynomial fit to the ^ statistic values near 
the minimum and determined the predicted minimum in 
order to improve the resolution of the fitting to f K, while 
saving computing time, a procedure which was verified 
by direct calculations. Figs. ^ and ^ show the best fits 
for the sun and Procyon from this method, including our 
final chosen spectral windows. Despite obtaining obser- 
vations of H7 of reasonable quality for some stars, this 
line is not considered in this work due to the problem of 
choosing suitable windows and the limited validity of the 
self-broadening calculations. 

In using such a method there is a natural tendency that 
a chosen window may not be entirely free from blends, 
as wings of nearby lines may encroach on the apparent 
window or unseen weak blends exist. One must therefore 
expect that the fitting method may very slightly over- 
estimate Toff in a systematic manner. This is likely of all 
such fitting, automated or not. 



4.1. Influence of spectral resolution 

It is important to understand the effect of spectral reso- 
lution on our results. Sufficient resolution is needed such 
that a suitable number of windows between blending lines 
are available. Once again, this is true not just of our fit- 
ting method but any approach. There must exist a limit- 
ing resolution for which all windows will be affected by the 
instrumental broadening of metal lines, but above which 
there is at least one unaffected window, and this limit dif- 
fers with star, and spectral region. For our purposes, a 
single small spectral window is insufficient, since as dis- 
cussed above, at the minimum, we would like to have sev- 
eral spectral windows evenly distributed across the Balmer 
line. 

Ideally, windows should be chosen such that for the 
test cases the derived effective temperature does not dif- 
fer between the high resolution spectra and the result ob- 
tained with the same spectra degraded to the resolution of 
our observations R « 30000, or at least introduces only a 
small error (say <10 K). This requirement was quite eas- 
ily fulfilled for Ha in the sun. However, for H/3 this was 
problematic. It proved difficult to find a set of windows 
which fulfilled both our criteria of a Teg in reasonable 
agreement with the high resolution result, and a reason- 
able coverage of the whole line profile shape. The windows 
eventually used for H/3 were chosen as a compromise be- 
tween both criteria, which resulted in a temperature from 
the degraded spectra of 17 K hotter than the result using 
high resolution spectra in both stars. This error is con- 
sidered acceptable, in view of gaining some information 
on the fit of the line shape, though we expect this error 
to increase in cooler stars. Because of this, and quite sig- 
nificant extra blending, H/3 will not be employed for the 
coolest stars in our sample. 



The results for different resolutions using the eventual 
set of spectral windows are incorporated in Table A 
brief survey of different resolutions indicates that a spec- 
tral resolution of at least R « 50000 is more appropriate 
for gaining both T^e and information on the line profile 
shape from H/3 for solar-type stars. 



5. Results 

5.1. Error Estimates 

Due to the complicated interplay of several mechanisms 
in the formation of hydrogen Balmer lines, estimating the 
sources of error and their behaviour with stellar parame- 



ters is similarly complicated. Fuhrmann et al. (1993, 1994) 
discuss many of these effects. Here, we have taken a quan- 
titative approach to the errors. By making a reasonable 
guess of the error in a number of possible sources, we 
have surveyed how these errors translate into errors in Teg 
through comparisons of model profiles. For consistency, 
the same windows as will be used for the observations are 
used, and test calculations employing the whole line pro- 
file indicate this does not introduce much error, which is 
also a confirmation that our windows do not introduce 
significant bias to certain parts of the line profile. Our 
error estimates are summarised in Table ^. Test calcula- 
tions showed that these errors do vary with MLT param- 
eters. For example, for the Tcff=6000 K and [Fe/H]=-2 
case. Table || indicates Ha to be quite sensitive to gravity 
while H/3 is rather insensitive. Test calculations for this 
case with MARCSa models in place of MARCSOSa find 
the sensitivity to gravity of Ha and H/3 to be roughly the 
same, a 0.1 dex change corresponding to about ±34 K for 
both lines. Therefore, our results are somewhat dependent 
on our choice of MLT parameters. 



Differences between the o ften u sed Vidal et al. ( 1973 , 
hereafter VCS) and Stehle ( |l994|) calculations perhaps 
give an indication of the magnitude of the error in the 
Stark broadening calculations, and these differences are 
reported in Table ^. Test calculations introducing a 5% er- 
ror to the Stark broadening found results of similar magni- 
tudes. VCS calculations always resulted in a stronger pro- 
file and therefore a cooler temperature. Stehle's (1994) cal- 
culations for proton perturbers are known to have reason- 
able agreement with experiment, better than VCS where 
ion dynamics are neglected. However, one should consider 
that while protons dominate in hot stars, in cooler stars 
the perturbing ions progressively become more a mix of 
protons and heavier ionised nuclei. Self-broadening for hy- 
drogen is less well studied than Stark broadening, no ex- 
perimental results currently existing. Based on error esti- 
mates for hydrogen broadening of metal lines using similar 
theoretical techniques we estimate the error at around 5% 



(e.g. Barklem & O'Mara 2001), though we note that the 
interaction is fundamentally different due to the resonance 
interaction. As the broadening by helium is only estimated 
here we adopt a rather large error bar of 50%. Typical er- 
rors in the gravity and metallicity are estimated at 0.1 
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Table 3. Derived Toff values for the sun and Procyon using different or degraded observations. Results are for 
MARCS05s models. Note all R and SNR values are approximate. 



Observation 


R 


SNR 


SNR 


T„ff Ha 


Toff RB 






Ha 


up 


fK) 


fK) 


Cliin 
kj UIl 












Kitt Peak FTS atlas 


> 3 X 10^ 


3000 


3000 


5733 


5723 




50000 


3000 


3000 


5735 


5729 




30000 


3000 


3000 


5739 


5740 


INT/MUSICOS (co-added) 


30000 


375 


240 


5743 


5748 


— (single exposure) 


30000 


170 


110 


5721 


5711 


Procyon 












McDonald/2dcoude 


2 X 10^ 




1000 




6474 




30000 




1000 




6481 


INT/MUSICOS (co-added) 


30000 


250 


160 


6538 


6487 


— (single exposure) 


30000 


120 


80 


6498 


6532 



dex. Possible interdependency of errors makes their cor- 
rect combination non-trivial, thus we chose a somewhat 
ad hoc procedure. The errors discussed to this point are 
grouped into two categories, broadening and stellar pa- 
rameters, and these two errors are combined in quadrature 
supposing they are independent and random. The totals 
for each category are found by simply summing. These 
fixed errors are sub-totalled in Table To estimate errors 
for a given line we interpolate in these fixed values which 
are then combined with the error sources discussed below 
which vary from case to case. 

As discussed in Sect. ^ observational errors are esti- 
mated at approximately 1%, and so the estimated errors 
in Table ^ corresponding to a 0.5% shift in the contin- 
uum placement represent perhaps a best case which we 
assume for our best spectra. In computing errors we in- 
terpolate in these values, and scale the result depending 
on the quality (SNR) of the observations employed. We 
assume 1.5% for our lowest quality spectra. In Sect. ^.2| 
we will see for our adopted model grid there is a discrep- 
ancy of order 50-60 K in the solar T^f value with the 
known value. Assuming negligible error in observations, 
gravity and metallicity, approximately half of this can be 
explained by our estimates of the error due to broadening 
theory given in Table ^. The remainder of the discrepancy 
can most likely be attributed to the models. Thus, for the 
solar case we estimate a conservative model error of or- 
der 40 K. It must be anticipated that this error will vary 
with stellar parameters, but it is difficult to estimate how. 
The error due to possible variation in MLT parameters is 
estimated from the difference between Tes values derived 
with MARCS and MARCS05 models, which is added to 
this 40 K. This above procedure amounts to using the 
sun to approximately calibrate our error bars. The effect 
of deviation from the solar or alpha-enhanced abundance 
patterns (i.e. abundance scatter) is estimated by half the 
difference between the results obtained from MARCS05s 
and MARCS05a model grids. We also include a fitting 
error, including resolution effects, of 20 K and 40 K for 



Ha and H/3 respectively, conservative values based on the 
tests we performed in Sect. 4.1 (see Table ||). Totals are 
shown in Table ^ for the case where there is no error due to 
convection or abundance scatter, and the best case obser- 
vations. Our final errors will account for these additional 
factors on a case by case basis. 



5.2. The sun 

The sun is the most important test case since high quality 
observations are available and the stellar parameters are 
well known. We follow Fuhrmann et al.'s ( |1993D approach 
and consider the MLT parameters as free and attempt to 
calibrate them at least approximately using the solar case. 
We say approximately since the appropriate convective ef- 
ficiency for representing Balmer lines must be expected 
to vary across the HR diagr am (e.g. Ludwig et al. 1999 ). 



Both Fuhrmann et al. (1993 ) using solar observations, and 
Steffen & Ludwig ( 1999 ) using hydrodynamical simula- 
tions of the sun, have found that a low value of a « 0.5 is 
needed to represent Balmer lines, both with y set at 0.5. 
We note that Steffen & Ludwig ( 1999 ) have shown using 
2D hydrodynamical models that, at least for flux-constant 
MLT models, "different mean [temperature] stratifications 
are needed to represent different spectroscopic properties 
of an inhomogeneous stellar photosphere" and therefore 
we should not expect the convective efficiency from obser- 
vations of Ha and H/3 which will be investigated here to be 
the same as that required to represent other observational 
quantities, even other Balmer lines. 

This work differs from Fuhrmann et al.'s ( 1993| ) in 
that, in addition to the mixing-length parameter a, the 
convective structure parameter y is allowed to vary, and 
the STEHLE+BPO broadening recipe is employed (see 
1993| for details of their broadening 



al. 



Fuhrmann et 

recipe) . We have investigated the solar case with the FTS 



solar atlas (Kurucz et al. 1984 ) using a number of mod- 
els where the MLT parameters a and y have been var- 
ied within reasonable limits. Two approaches are taken. 
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Table 4. Estimated errors in T^b (in K) corresponding to introduced errors in the input physics, stehar parameters 
or observations for MARCS05 models, with MARCS05a models used for [Fe/H]= —1 and —2 cases. Results are given 
for various Toff at both solar and low metallicity (quoted as Toff/pe/H]). All cases use \ogg = 4.2, except the sun 
where logg — 4.44. Errors are introduced (amounts discussed in the text) and the corresponding changes in derived 
Teff observed. Signs show the direction of the change, for the corresponding change in the input physics or stellar 
parameter. Errors are usually symmetric, but if not the largest is quoted. We emphasise these errors vary with chosen 
MLT parameters. For a given stellar spectrum the continuum error is interpolated from the listed values and scaled 
depending on the SNR, and the values of AToff (convection) and ATcff (abundance scatter) are computed case by case 
and given in Table ll 



Introduced 


Sun 


5000/0.0 


6000/0.0 


7000/0.0 


5000/-1.0 


6000/-1.0 


7000/-1.0 


5000/-2.0 


6000/-2.0 


7000/-2.0 


error 


Hq/H/3 


Hq/H/3 


Ha/H/3 


Hq/H/3 


Ha/H/3 


Hq/H/3 


Hq/H/3 


Hq/H/3 


Hq/H/3 


Hq/H/3 


Stark-broadening 






















VCS 


-5/-15 


-3/-10 


-7/-18 


-19/-28 


-1/-7 


-7/-17 


-18/-26 


0/-3 


-8/-17 


-17/-25 


Self-broadening 






















±5% 


T14/T7 


T10/T6 


T13/T6 


T8/T3 


T15/T13 


T23/T11 


T11/T3 


=F23/=f20 


=F38/t15 


T12/T4 


Hc-broadcning 






















±50% 


T6 /t4 


T5/T3 


T6/T3 


T4/T2 


T6/T7 


T10/T6 


T5/T2 


TlO/Tll 


T17/T8 


T6/T2 


logs 
±0.1 dex 




T6/T6 


T9/T3 


T3/±8 


T13/T12 


T25/T7 


T3/±ll 


T33/T17 


T52/T11 


T5/±ll 


[Fe/H] 






















±0.1 dox 




±19/-2 


±7/Tl7 


=F25/=f31 


±35/±20 


±17/=f7 


T14/T15 


±16/t7 


±6/t4 


±6/±6 


Sub-total 


±25/±26 


±31/±21 


±31/±34 


±42/±51 


±53/±42 


±58/±37 


±38/±40 


±59/±48 


±86/±43 


±37/±35 


Continuum 






















±0.5% 


±38/±29 


±36/±23 


±29/±19 


±34/±21 


±33/±25 


±38/±23 


±35/±21 


±38/±25 


±56/±26 


±36/±21 


Model 










- 40 + |AT„ff (convection 










Abundance scatter 










ATcff (abundance scatter) 1 > 








Fitting &: resolution 












20/40 > 










Total (best case) 


±69/±72 


±65/±65 


±61/±69 


±70/±79 


±76/±75 


±82/±71 


±69/±73 


±83/±79 


±lll/±76 


±68/±70 



Firstly, we derive Toff considering it a free parameter for 
a given MLT a, y set and record the value. Secondly, 
we set the effective temperature at 5777 K and survey the 
value across the different MLT parameters. The results 
are displayed in Fig. ^, where all values are normalised 
so the best possible fit for a given line has = 1- Steps 
in of 0.5, from 1 to 20, arc plotted (except in one case 
where we use steps of 1 for clarity) . Differences below this 
level correspond to only small differences in the line pro- 
files and are probably meaningless. The contour at 1.1, 
where it exists, is also shown to illustrate the plateau at 

It is seen that no single parameter set shows itself 
to be especially preferred on the basis of agreement with 
line shape and Teg. A simple approximate picture is that 
agreement with the line shape indicates correct temper- 
ature structure in the line formation region, and agree- 
ment with Teff indicates correct temperature at the line 
formation region, although the reality is somewhat more 
complex. The line profile shapes for both Ha and H/3 
are best fit anywhere on the plateau at low a and low 
y with a T^s some 50 K lower than the accepted value. If 
we force Toff = 5777 K the fits to line profile shapes for 
both Ha and H/3 are quite poor, typically being too deep. 
Furthermore, a single a, y combination will not give the 
correct T^s for both lines, seen by the fact that the loci 
of a, y values giving the correct solution (thick lines) do 
not overlap. More generally one can see from Fig. || that 
Balmer lines do not well constrain the MLT parameters. 



since an increase in a can be compensated by a decrease 
in y, further strengthening Steffen & Ludwig's (1999) as- 
sertion that Balmer lines should not be seen as evidence 
for low efficiency of solar-type convection. 

On the assumption that the FTS observations em- 
ployed are of high accuracy in comparison with the 
model spectra, one interprets this as a shortcoming of 
either the model atmospheres or the broadening the- 
ory, though most likely a combination of both. In pa- 
per II using the same broadening theory we found that 
the semi-empirical model of Holweger & Miillcr ( 1974 ), 
which matches limb- darken ing better than MARCS mod- 
els (Blackwell et al. 1995| ), in fact gives lines too shal- 
low. Therefore it is certainly feasible that this discrepancy 
could be due to models. Based on these results we decided 



to adopt the same parameters as Fuhrmann et al. (1993) 
namely a=0.5 and y=0.5 as our default values since this 
parameter combination reproduces the line shapes and 
therefore the temperature stratification in the line form- 
ing region as well as any other choice. As we see in Fig. ^ 
the fit to H/3 is not perfect showing differences of order 
1% with a trend across the line, again indicating some de- 
ficiency in either models or line broadening. Although a 
lower value of y is perhaps more justified physically on the 



basis of the diffusion approximation (Henyey et al. 1965), 
adopting t hese param eters will make comparisons with 
Fuhrmann (|199^ , ^OOOl) more straightforward. In any case, 
all calculations will be repeated with a=1.5 and j/=0.076 
as a test of sensitivity to this choice. The error in T^ft in- 
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troduced by this choice of MLT parameters was discussed 
in Sect. 



5.1 



In addition to testing models, the sun is also the best 
test case for our observational data and reduction proce- 
dures. By applying our fitting method to different obser- 
vations we can test their accuracy. A high SNR spectrum 
of the moon was produced from co-added spectra, as well 
a spectrum at lower SNR more typical of our other ob- 
servations from a single exposure. The two spectra were 
reduced completely independently. Comparison of the de- 
rived Teff values with that from the FTS atlas degraded to 
the same resolution presented in Table ^ shows a scatter 
of typically ±15 K for Ha and ±30 K for H/3. 

5.3. Procyon 

The F5 star Procyon (HR 2943) is of fundamental impor- 
tance as it is one of the few dwarf stars other than the sun 
accessible to key direct measurements due both torts prox- 
imity and binary nature. AUende Prieto et al. ( 2002 ) find 
Tcff = 6530 ± 49 K and log.g = 3.96 ± 0. 02 in g ood agree- 
ment with other works (e.g. Di Benedetto 1998, Fuhrmann 
et al. [1997[ ). 

As mentioned, we have a high quality spectrum of H/3 
for Procyon, and therefore this star serves as a second 
important test case. Results shown in Table |^ indicate 
that for our window set R ~ 30000 is adequate for H/3, a 
result which is expected since blending is less pervasive. 
When comparing spectra from different observations at 
R 30000, a scatter of about ±50 K is seen. This is more 
than for the solar case, which we interpret as being due 
to the lower SNR of the single exposure leading to a less 
reliable continuum placement. However, we should also 
consider the possibility that the broader lines in Procyon 
make the reduction procedure slightly less intrinsically re- 
liable. 

For this star, with the best available spectrum, we find 
6538±82 K from Ha and 6474±143 K from H/3. Our final 
adopted result (weighted average, rounded) of 6520±80 K 
is in good agreement with the above mentioned results. 
We note also that this result agrees within error with 
IRFM results of 6579±100 K from Alonso et al. ( |1996[ ) and 
6601±75 K from Saxner & Hammerback ( 1985 ), bu t does 
not agree with the value found by Edvardsson et al. (1993) 
of 6704 K from h — y colours. We note, that as for the solar 
case the fit to H/3 shown in Fig. || is not perfect showing 
differences with a trend across the line profile. 



5.4. The sample 

Using our automated fitting procedure, effective tempera- 
tures were derived from all four model grids (MARCSs, 
MARCSa, MARCS05S, MARCS05a) using both broad- 
ening recipes, adopting surface gravities and metallicities 
from the literature. Fitting for all four models with both 
broadening recipes allows some interesting comparisons, 
and is also useful in estimating errors as discussed in 




Fig. 7. The difference between derived Teff from Ha and 
H/3 using MARCS05 and STEHLE+BPO. 



Sect. 5.1 . Following the well known enhancement of alpha- 
elements, for stars with metallicity less than [Fe/H]= —0.5 
the result of the alpha-element enhanced grid is used, oth- 
erwise the scaled-solar result is adopted. Unless otherwise 
stated, all discussion below refers to the STEHLE+BPO 
results. 

The derived values of Teff from MARCSG5 for each 
line are presented in Table ^ the adopted temperature 
being the weighted average of the results from both lines. 
Some selected example fits to the spectra are shown in 
Fig. ^. Fig. plots the difference between T^s from Ha 
and H/3 against metallicity (no trend is seen with Teff), 
and suggests that H/3 gives temperatures which are cooler 
by 69±70 K than those from Ha. Essentially the same 
result is found for STEHLE+AG recipe. When MARCS 
models are employed we find better agreement in that H/3 
is typically cooler by 11±75 K, though the line shapes 
become noticeably worse. Comparing the fits to H/3 for 
MARCS and MARCS05 models it was found that on av- 
erage Xmarcs/Xmarcsos = 1.85±0.75, all stars indicating 
a better fit to the shape with MARCSG5. Ha also showed 
better fits for MARCS05 for most stars although the ef- 
fect was not as strong, as would be expected due to the 
decreased sensitivity of Ha to the convection parameters. 
Very similar results were obtained with the STEHLE+AG 
recipe. One motivation for using a quantitative profile 
comparison was to look for trends with stellar parameters 
in the relative goodness-of-fit of model profile shapes from 
different convection parameters to observations. However, 
we note that the Xmarcs/Xmarcsos ratio is expected to 
be dependent on SNR, which was confirmed from experi- 
ments using solar and Procyon spectra. Although the re- 
sult converges for high SNR, at the typical SNR used here 
it was found the ratio would be significantly underesti- 
mated. Therefore care must be exercised in comparisons 
of such statistics between different stars and models. 
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Table 5. The programme stars and derived effective temperatures. Sources indicates tlie source of adopted logg and 
[Fe/H] values in that order. Estimated errors due to convection and abundance scatter are shown, where ATeff(conv)= 
Toff(MARCS)-roff(MARCS05) and ATeff(abund)= ±0.5 x (reff(MARCS05a)-Teff(MARCS05s)) with the negative 
case of this equation corresponding to [Fe/H]< —0.5. Signs are shown to indicate the direction of the change, but the 
error is the absolute value. The error bars shown for the sun are if it is treated as all other stars, despite that a few 
of the errors are much less for this case. 



Star 


logs [Fe/H] 


Sources 


Tcff 




ATeff 




ATeff 


ATeff 


Teff 


HR HD 


adopted 




Ha 


Ha 


Ha 


H/3 


H/3 


H/3 


adopted 










conv 


abund 




conv 


abund 






(cgs) 




(K) 


(K) 


(K) 


(K) 


(K) 


(K) 


(K) 



Sun 




4.44 


0.00 




5733± 81 


-11 


39 


5723± 90 


19 


9 


5730± 70 


321 


6582 


4.39 


-0.86 


AFG, AFG 


5449±100 


6 


-58 


5288±135 


57 


-24 


5390± 90 


753 


16160 


4.50 


-0.08 


A AM, A AM 


4843±124 


3 


52 








4840±130 




19445 


4.50 


-1.88 


AP, GCC 


5783±172 


-29 


-18 


5849±121 


33 


-1 


5830±100 


996 


20630 


4.50 


0.12 


BLG, BLG 


5733±103 


-9 


38 


5682±102 


19 


7 


5710± 80 


1101 


22484 


4.02 


-0.11 


AP, GCC 


6033± 79 


-16 


24 


6008±107 


40 


-3 


6020± 70 




22879 


4.38 


-0.76 


AP, GCC 


5867±119 


-5 


-45 


5765±126 


45 


-13 


5820± 90 


1325 


26965 


4.65 


-0.21 


BG, BG 


5187±111 


1 


57 








5190±120 


2852 


58946 


4.17 


-0.31 


SH, SH 


6959± 73 


-5 


-6 


7006±184 


127 


-12 


6970± 70 


2943 


61421 


3.96 


-0.06 


APP, AP 


6538± 82 


-22 


4 


6474±143 


82 


-13 


6520± 80 




64090 


4.66 


-1.60 


AP, GCC 


5342±143 


4 


-54 


5175±166 


85 


-24 


5270±110 


3775 


82328 


3.61 


-0.01 


AP, AAM 


6302±101 


-23 


9 


6205±117 


53 


-6 


6260± 80 


4150 


91752 


3.84 


-0.21 


GCC, GCC 


6567±116 


-24 


1 


6404±135 


69 


-9 


6500± 90 


4540 


102870 


3.98 


0.10 


AP, GCC 


6140±100 


-15 


19 


6021± 98 


26 


-1 


6080± 70 


4550 


103095 


4.68 


-1.22 


AP, GCC 


5113±108 


17 


-59 








5110±110 


4623 


105452 


4.00 


-0.60 


AAM, AAM 


7029± 97 


-11 


6 


6996±156 


96 


8 


7020± 90 




108177 


4.47 


-1.55 


AP, GCC 


6071±171 


39 


-11 


6087±148 


75 


3 


6080±120 


4983 


114710 


4.32 


0.00 


AP, GCC 


6025± 80 


-15 


32 


5971±101 


37 


2 


6000± 70 




114762 


4.22 
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Sources for log q and [Fe/H] parameters: A AM — 
Prieto et al. ([l999[), AP P = A Uende Prieto et al. (|oo|) 



Alonso et al. ( |l99q ), AFG = Axer et al. (|1994|), AP = AUend e 
BG = Bellfc Gustafsson (|l98£|), BLG = Blackwell & Lynas-Gray (|l994|). 



GCC = Gratton et al. (1996), SH = Saxner & Hammarback (1985) 



6. Comparisons and discussion 

Fig. compares our adopted effective temperatures 



against the results from the IRFM of Alonso et al. (1996) 
for both STEHLE+AG and STEHLE +BPO recipes. We 
choose comparison with Alonso et al. ( |1996D as their re- 
sults include both solar metallicity and metal-poor stars, 
and has the largest overlap with our sample. Using 
the STEHLE-I-AG recipe we find that the tempera- 
tures are on average hotter than those found from the 
IRFM by 110±77 K. When the STEHLE+BPO recipe 
is used we see generally agreement within error, the 



results typically hotter by 11±95 K. Two outliers are 
seen, namely HD 19445 and HD 64090, which we find 
to be signific antly cooler than the IRFM values from 
Alonso et al. ( 1996 ). If we compare only those 4 stars in 
our sample with metallicity below —1.5, which includes 



these two stars, with Alonso et al. (1996) a difference 
of -117±97 K is found. The only other IRFM calibra- 
tion we are aware of in the literature for metal-poor stars 
is by Magain ( |1987| ), who finds temperatures typically 
112±56 K cooler than Alonso et al. (1996) for 11 metal- 
poor stars. Alonso et al. ( |1996D attribute this zero point 
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difference to tfie atmospheric models used. If we again con- 
sider only the 4 stars in our sample with metallicity below 



— 1.5 and compare with the results of Magain (1987) the 
difference is — 34±85 K, our temperatures from Balmer 
lines being slightly cooler. These 4 stars are the only com- 
mon stars with the Magain (1987) sample. Thus, our re- 
sults are in better agreement with Magain (1987) for the 
metal-poor stars, though 4 stars is insufficient for any clear 
conclusion. 

There is a suggestion from Fig. ^that our t empe ratures 
are systematically hotter than Alonso et al. ( 1996 ) in the 
coolest stars in the sample. This may well be due to the 
increased blending causing an overestimate of the Balmer 
line temperatures, though we note the IRFM may also be 
less reliable in this regime. 

Similar comparisons were made for common stars with 
the results for population I stars of Blackwell fc Lynas- 



Gray (|l99j, |l998| ) and Saxner & Hammerback (|l985| ) both 
using the IRFM and the results showed generally agree- 
ment within error. Agreement with T^s values for 16 com- 
mon stars from Edvardsson et al. (1993) using b — y colours 
was also good, with a mean difference of —51 ± 67 K. 
Two exceptions were Procyon and HR 2852 both early F 
stars, which we found to be significantly cooler than the 
Edvardsson et al. (1993) temperatures. 

Fuhrmann et al. (1994) determined temperatures for 
a large sample of stars using Balmer lines, these de- 
terminations having been superseded by more recent 
work with better observations and improved models 
(Fuhrmann 1998, 2000). We compare our temperatures 



with Fuhrmann's work (1998, pOOO ) in Fi g. P (r ight panel). 
As Fuhrmann employs the Ali & Griem (1966) theory we 



might have expected to see differences of the order pre- 
dicted in paper II (Fig. 7 of that paper). In fact, the agree- 
ment around solar metallicity and even slightly below is 
good. Comparison of our results using the STEHLE-fAG 
recipe with those of Fuhrmann are shown in Fig. I (left 
panel). Reasonable agreement between the two sets of 
results might have been expected, but the results show 
a systematic difference of 91±40 K with no significant 
trend with stellar parameters. The reason for this dif- 
ference has been traced to a number of factors, namely 
systematic differences in Stark broadening, model atmo- 
spheres, and a slightly different application of the Ali & 
Griem (1966) theory (see Sect. |3|). These differences are 
all quite small but by chance all act in the same direction, 
with relative importance varying with stellar parameters 
and line. The differences are also blurred by random errors 
such as the use of different observations, fitting methods 
and adopted gravity and metallicity. For example using 
MARCS05 and STEHLE-hAG we fit the solar Ha s pectra 
best with Toff « 5830 K, while Fuhrmann et al. ( [1997 ) 
fit the line with T^ff w 5750 K. As shown in Sect. 3.1, 
« 20 K of this originates in the atmospheric models. 
Comparisons of the Stark broadening tables employed by 
Fuhrmann et al. were found to give very slightly stronger 
profiles than VCS, and thus an even greater difference 
when compared to STEELE than VCS. This accounts for 



a further « 10 K (Table If one adds the resonance 
broadening of the 3p state to that of the 2p state, the 
broadening is increased by about 15% compared to con- 
sidering only the 2p state as we do in our STEHLE-fAG 
recipe. We see from Table |^ that this produces a signifi- 
cant difference of « 40 K. As stated above all these dif- 
ferences act in the same sense in that Fuhrmann et al.'s 
profiles are stronger, and thus explain the difference of 
« 80 K in derived temperatures. Similarly, for Hf] we ob- 
tain Toff ~ 5820 K, and Fuhrmann et al. Toff « 5750 K. 
For this line the difference in models and implementation 
of resonance broadening is less important, w 10 K each, 
while the Stark broadening difference is increased ~ 20 K, 
explaining w 40 K of the difference. The remainder is cer- 
tainly within fitting and observational error. We note that 
these differences of order 80 K between the two determi- 
nations are consistent with, and thus support, our error 
estimates. 

7. Conclusions 

Effective temperatures have been derived from Balmer line 
profiles using a quantitative fitting method with a detailed 
error analysis including investigation of the susceptibility 
to various errors with temperature and metallicity for our 
models. Our temperatures find good agreement with the 
IRFM near solar metallicity but show differences at low 
metallicity where the two available IRFM determinations 
themselves are in disagreement. Our results for metal- 



poor stars seem in better agreement with Magain (1987), 
though we caution this is only for 4 stars and the better 
agreement may be fortuitous. This should be investigated 
further, for which more, and better observations of Balmer 
lines in population II stars are needed. The origin of the 
differences between these two IRFM determinations needs 
also to be understood. 

Through estimates of the uncertainties, we found that 
the relative weight that should be given to Ha and II/3 in 
determining effective temperatures varies quite substan- 
tially with stellar parameters. This is predominantly a 
matter of balance between the reduced sensitivity of Ha 
to temperature at low metallicities making broadening, 
logg and observational uncertainties very important, and 
the high sensitivity of H/3 to convection. Thus to improve 
the accuracy of Balmer line temperatures for metal-poor 
stars these four areas must be addressed. Observational 
uncertainties will continue to improve as large telescopes 
provide the possibility of combined high resolution and 
SNR for these faint objects. There no doubt exists a lim- 
iting accuracy for continuum determination procedures for 
echelle spectra such as that used here, but we believe this 
has not yet been reached. Fuhrmann et al. (1994) sug- 



gest an accuracy of about 0.3% is achievable with mod- 
ern spectroscopy. Planned satellite astrometry missions 
should address the gravities. Both the remaining uncer- 
tainties, convection and self-broadening, lie in the theo- 
retical realm at least at present. Considerable progress 
has been made in 3D hydrodynamical simulations of con- 
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vection and inhomogeneities in cool stars, and this will be 
the way forward. However, for the immediate future where 
direct application of such models is impractical, calibra- 
tions of MLT parameters for Balmer lines across the HR 
diagram similar to Ludwig et al. ( |1999D or tabulated Tcff 
corrections (against a given MLT parameter set) would be 
important. Improved calculations of the self-broadening of 
Balmer lines should be undertaken without resort to the 
impact approximation which will require improved short 
range potentials from those used in paper II. 

In paper I and II large differences between Toff val- 
ues derived with Ali & Griem (1966) theory and new 
calculations of the self-broadening, particularly in metal- 
poor stars, were predicted. Comparison of our results with 
those of Fuhrmann ( |1998| , |2000D where Ali & Griem the- 
ory has been employed, find in fact good agreement ex- 
cept for low metallicity stars of around solar tempera- 
ture where the differences are smaller than expected, but 
the temperatures are still significantly cooler. The reason 
for this unexpected agreement is a number of systematic 
differences in models and employed broadening recipes, 
which though typically individually small, together some- 
what compensate the difference made by the new self- 
broadening. This emphasises that in order to achieve high 
absolute precision in T^s determinations from Balmer lines 
such small differences must in fact be carefully considered. 
Errors in the relative temperature scale, particularly be- 
tween stars of low and solar metallicity important in trac- 
ing chemical evolution in the galaxy, are unlikely to be 
much smaller than the errors in the absolute scale since 
those errors which might cancel in differential compari- 
son, namely broadening theory and model errors (includ- 
ing convection), typically vary with metallicity resulting 
in little cancellation. 
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Fig. 5. Survey of Balmer line fits for the observed solar fliix spectra with different MLT parameters using 
STEHLE+BPO recipe. The left column of plots refers to Ha and the right to H/3. In each case the uppermost 
plot is the values for the best fit in the case where Tgs is considered a free parameter, the thick line on these plots 
indicating the locus where Toff = 5777 K is found. The middle plot shows the corresponding behaviour of the best fit 
Tgff values. The lowermost plots show for the case where Tcs is fixed at 5777 K. All are normalised so the best 
fit for a given line is 1.0, (i.e. both top and bottom plots of each column are on the same scale) which in both cases 
corresponds to the plateau in the uppermost plots at low a and y. Filled circles show the points actually computed 
and the contours are found by interpolation (minimum curvature surface) and thus any structure on scales smaller 
than the computed grid should be ignored. 
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Fig. 6. Selected Balmer line spectra and fits. Full lines are the observed spectra and dotted lines the synthetic spectra 
corresponding to our best fit with MARCS05 and STEHLE+BPO, the parameters for which arc given to the upper 
right of each profile (Teg / log <?/ [Fe/H] ) . Tick marks on the vertical scale correspond to . 1 in residual flux, and continua 
are at every 3rd tick mark. Double-thickness dashed vertical lines show the approximate validity region of the impact 
approximation of the BPO theory (for Ha this always lies off the plot and is therefore only seen for H/3). Shaded 
regions show windows used for fitting for each case. 
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Fig. 8. Derived temper ature s for STEHLE+AG (left) and STEHLE+BPO (right) recipes compared with the IRFM 
results of Alonso et al.( 1996| ). 
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Fig. 9. Derived temperatures for STEHLE+AG (left) and STEHLE+BPO (right) recipes compared with the Balmer 
line results of Fuhrmann against metallicity. 



